Spin-1/2 geometric phase driven by decohering quantum fields 
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We calculate the geometric phase of a spin-1/2 system driven by a one and two mode quantum 
field subject to decoherence. Using the quantum jump approach, we show that the corrections to 
the phase in the no-jump trajectory are different when considering an adiabatic and non-adiabatic 
evolution. We discuss the implications of our results from both the fundamental as well as quantum 
computational perspective. 
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In quantum mechanics physical states are equivalent 
up to a global phase which in general does not contain 
useful information about the described system, and thus, 
can be ignored. However, Berry surprisingly showed 
that these phases can have a component of geometric ori- 
gin with important observable consequences, being the 
most cited examples the Aharonov-Bohm effect |2| and 
the spin-1/2 particle driven by a rotating magnetic field 
0. These components which are gauge invariant and 
only depend on the path followed by the system during its 
evolution, have been investigated and tested in a variety 
of settings and have been generalized in several directions 
[j| . Geometric phases are interesting both from a funda- 
mental point of view and for their applications, among 
which geometric quantum computation [J, |3, |6j is one of 
the most important. In fact, the use of geometric phases 
in the implementation of fault-tolerant quantum gates 
has motivated their study under more realistic situations 
0, 0- For example, when a system interacts with an 
environment, its quantum superpositions may decay into 
statistical mixtures (SlllfJ and this effect, called decoher- 
ence, is the most important limiting factor for quantum 
computation. 

Previous works investigate the behavior of geomet- 
ric phases under some typical errors sources like ran- 
dom classical fluctuations to the driving fields, as well as 
generic reservoirs acting in spin 1/2 evolutions 0, • 
All of them consider the driving field as a classical sys- 
tem. However, any driving field is also a quantized sys- 
tem and, whenever this quantum behavior is relevant, 
which is the case in many experimental situations, deco- 
herence on these fields may play an important role. In 
fact, it may even become critical, particularly when geo- 
metric phases are used to implement quantum protocols, 
like communication and computational ones. 

In this letter, we investigate the behavior of the geo- 
metric phase of a spin 1 /2 particle interacting with a driv- 
ing magnetic field when this field is not only quantized 
but also subjected to decoherence. We calculate and ana- 
lyze the effect of decoherence of the driving field on both 
adiabatic and non-adiabatic evolutions of the joint spin 
and quantized modes system. First we briefly describe 



the general framework of geometric phases in open sys- 
tems, developed in Then we calculate Berry's phases 
for different interactions of spin 1/2 systems and decoher- 
ing fields both in the adiabatic and non-adiabatic scenar- 
ios. Finally we point out the differences between these 
two situations and how this noise source compares to pre- 
viously analyzed ones. 

Let us first consider a system described by the den- 
sity operator p and a Hamiltonian H. The decoherence 
process due to the interaction with an environment (un- 
der the Markovian approximation) is described by the 
following master equation (h = 1): 
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where the commutator generates the coherent part of the 
evolution and the remaining part represents the effect of 
the reservoir on the dynamics of the system. The action 
of each amounts to a different decohering process. 
Suppose that we monitor the system and do not detect 
any decay. The geometric phase for the "no-jump" tra- 
jectory for master equation in the continuous limit is 
given by |l2j |: 
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where 



i-\<p°(t)) = H\r(t)), hA») = hM, (3) 

and H is a non-Hermitian effective Hamiltonian given by: 
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Before applying this general framework to our problem 
we consider the adiabatic phase of a two-level system and 
a single mode quantum field following We describe 
the two-level system with Bohr frequency u) in terms of 
Pauli operators a z , cr± = (cr x ± *Oj,)/2 and the field with 
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frequency v in terms of the creation and annihilation op- 
erators a and a'. In the interaction picture, the initial 
Hamiltonian reads 



Hint = ycr 2 + g{da- + acr + ) 



(5) 



where A = lu — v is the detuning between the quan- 
tum mode and the two-level system and g is the coupling 
constant. The evolution of the system is dictated by 
the usual time dependent Schrodinger equation. Follow- 
ing 0|, to generate a geometric phase the Hamiltonian 
is varied in a cyclic and adiabatic fashion by means of 
the unitary operation U = e~ l ^ n where h — a^a is the 
number of photons in the field. Using Berry's formula 
for the phase Q we find that the eigenstate of © 

|*+) =cos(0„/2)|e,n)+sm(0 n /2)| 5 ,7i+l), (6) 



acquires the phase 



7 + = 2n{n) = 2nn + tt(1 — cos 9 n ), 
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where cos#„ = ^ with R = -\/A 2 /4 + g 2 (n + 1), when 
<p is varied from to 2ti. To consider now that the field is 
subject to decoherence we describe our system by equa- 
tion Q with r = ^/Xa, i.e. the field is linearly losing 
photons to its reservoir. Thus, the effective Hamiltonian 
is H = H — i^h. Using equation and © we obtain 
the phase 
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This result can be easily derived using the bi-orthogonal 
basis technique illustrated in [J^ which gives rise to a 
complex geometric phase whose real part corresponds, in 
the adiabatic case, to the general formula 

Notice that the lowest order correction in X/R of eq. (JSJ) 
is quadratic, which means that in case of low decoherence 
we recover the phase 7+, up to the first order in X/R. To 
the second order in X/R the geometric phase (JSJ reads: 



7+ + 7T cos ( 
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(1 + 3/2 cos 2 9 n ) . 



It is important to notice that in case of low decoherence 
the deviation of the geometric phase from the value 7+ 
is null up to the first order in X/R. This reflects the 
resilience of the geometric phase against the environment. 
The reason for this can be interpreted heuristically as 
a consequence of the adiabatic evolution of the system 
under consideration. 

Under the adiabatic assumption, the state of a sys- 
tem, in its evolution, tends to follow the eigenstates of 
the instantaneous Hamiltonian. Loosely speaking, this 
behaviour opposes to the tendency of the environment of 
dragging the state away from its undisturbed evolution. 



More precisely, the adiabatic approximation ensures 
that the probability for a state to follow the instanta- 
neous eigenstate of the unperturbed Hermitian Hamilto- 
nian is highly enhanced compared to the probabilities of 
transitions to other eigcnspaces. In fact, the probability 
amplitudes associated to these transitions are averaged 
down due to their high frequency evolutions (of the order 
of Bohr frequencies of the system) , whereas the probabil- 
ity amplitude of staying in the same eigenspace is almost 
stationary. If the decoherence rate is sufficiently small, 
this amounts to effectively projecting the state back to 
the original eigenspace, whenever it tends to be driven 
away by the environment. Therefore, the state trajectory 
on the projective Hilbert space tends to be unaffected by 
the decoherence (up to the first order), thereby leaving 
the area enclosed in the path, and hence the geometric 
phase, unchanged. 

This characteristic of robustness against decoherence 
is also present in other systems. For example, similar 
properties can be observed in the analogous model, an- 
alyzed in ref. |l3|. in which a two-level atom interacts 
with not one but two quantized modes of a harmonic os- 
cillator. In this case, the geometric phase is obtained by 
an adiabatic evolution of the initial Hamiltonian 



H = va' a + vb^b + g(<7+a + cr_a^), 



(9) 



where the extra mode is described by the creation and 
annihilation operators b and respectively. 

Analogously to the previous case, an initial eigen- 
state of this Hamiltonian \ip) — cos #„/2|e, n, n') + 
s'm9 n /2\g,n + l,n'), is considered. After a cyclic evo- 
lution of this state following the adiabatic rotation of 
Hamiltonian JHJ in a two-dimensional parameter space, 
the original state acquires a geometric phase equal to: 



X(n,n>) = 2° 



n — n' + i (1 — cos n ) 
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where fi is the solid angle described by the parameters. 
In ref. [l5| the authors suggest a cavity QED experiment 
that implements this rotation and allows for the mea- 
surement of the above mentioned phase. 

In order to study the case where the interacting fields 
are decohering, we consider the Hamiltonian H = H — 
i^N with N — a) a + b^b the total number of photons 
in the system. As in the single-mode case, here a non- 
Hermitiam Hamiltonian is obtained from the assumption 
that no jump occurs during the evolution, i.e. the system 
is assumed to be continuously monitored by detectors and 
no emission of photon is registered. The adiabatic phase 
then yields 



(11) 
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Again, the expected geometric phase (fTT))) is recovered in 
case of low decoherence, and the lowest order correction 
is only quadratic in X/R, being R = ^/A 2 /4 + g 2 (n + 1) 
the Rabi frequency. Finally, note that a second order 
correction in the decaying factor A for the "no-jump" tra- 
jectory suggests that the fields decoherence may not play 
such an important role in the realization of the proposed 
experiment 

It is interesting to compare the geometric phase due 
to an adiabatic evolution in presence of decoherence and 
the analogous result obtained in a non-adiabatic fashion. 
It is well known ^(| that the adiabaticity is not a nec- 
essary condition to observe geometric phases. In fact, 
these are uniquely defined by the path on the projective 
Hilbert space traversed by the quantum system in its evo- 
lution. Thus, no matter how this evolution is achieved, 
the geometric phase will remained unchanged. This may 
no longer be the case in presence of decoherence, as the 
interaction with the environment can affects differently 
adiabatic and non-adiabatic evolutions. 

A typical example of non-adiabatic evolution is the one 
in which the Hamiltonian is time-independent and the 
initial state is chosen to be a cyclic state, i.e. a state 
that after a suitable time T evolves back to itself, up to 
a phase change. In some cases it is possible to choose 
a Hamiltonian and a cyclic state such that they gener- 
ates exactly the same path as the one associated to a 
given adiabatic evolution. In |l7| we considered the non- 
adiabatic version of the fully quantised spin-1/2 phase. 
Here we present the analysis of this case taking into ac- 
count the effect of decoherence and compare the geomet- 
ric phase obtained in the analogous adiabatic case. 

In the non-adiabatic setup that we are presenting, the 
system is initially prepared in an entangled state of atom 
and the field mode. Then it evolves under a time indepen- 
dent Hamiltonian involving only the degrees of freedom 
of the field. By turning on a Jaynes-Cummings interac- 
tion we can prepare the system in the state (0- 

After this interaction has been switched off, we assume 
that the dynamics of the system is described, in the in- 
teraction picture, by the Hamiltonian: 

H int = pn, (12) 

where (3 is a constant parameter. Thus, the state evolves 
according to: 

|*+(t)) = e-^ ftt |# + ), (13) 

and after a time T = 4ir/f3 the state completes a closed 
loop. Using the definition of Aharanov and Anandan 
geometric phase 0] it is easy to show that after a cyclic 
evolution the phase acquired by the state is 7+, the same 
as © obtained in the adiabatic case. 

When the decoherence of the field is considered, the 
phase for the no-jump trajectory can be calculated from 



the expression J5J) which yields 

/ = -^ln(^+| e - 27r W fi |*+). (14) 
A 

Given the initial state (JBJ, the eq. (|14|1 reads 

7^ = 2nn - f ln(cos 2 — + e~ 2 ^ sin 2 — ), (15) 
A 2 2 

As expected, the decoherence-free case is recovered for 
low values of the parameter A. In fact in the limit of 
A << (3 the geometric phase results: 

l c l = 7 + + ^(™#„) 2 + ((^) 2 ). (16) 

This result can be easily understood in the following way. 
In the case of no-jump evolution, i.e. of no photon emis- 
sion, the decay rate due to the imaginary part of the 
Hamiltonian H is proportional to the number n of pho- 
tons contained in the mode. Therefore, given an initial 
superposition of two states with different values of n, 
the amplitude associated to the lowest number of pho- 
tons gradually increases in time and eventually the sys- 
tem converges to the state with the lowest n. In other 
words, since, for the no-jump trajectory, no photon de- 
caying is observed, the probability of the lower n state 
increases with time. In the Bloch sphere representation 
of the subspace {|e, n), \g, n+1)}, the evolution under the 
free Hamiltonian i|18|) would just appear as a rotation of 
the Bloch vector around the z axis. By considering the 
decoherence of the field the state will then spiral towards 
the south pole. Thus, the first order term in A/ (3 appear- 
ing in Eq. (|16|l accounts for the extra area spanned by 
the system on the Bloch sphere due to the decoherence. 

As expected the geometric phase is affected by decoher- 
ence in different ways for the adiabatic and non-adiabatic 
scenarios. In particular, for low decoherence rates, i.e. 
A << R and A << (3 in the adiabatic and non-adiabatic 
case, respectively, the lowest correction is quadratic in 
the former and linear in the latter case. This should be 
expected since in the adiabatic evolution the probabil- 
ity for the state to be dragged away by the decoherence 
from the unperturbed evolution is washed away (in the 
first order) by the driving Hamiltonian, thereby oppos- 
ing against decohering effects. On the other hand, in the 
non-adiabatic evolution, there is no action other than the 
decoherence, which finds no resistance in the evolution. 

Analogous considerations can be done in the case of 
a two modes system described by the Hamiltonian 
The non-adiabatic version of the same problem can be 
described in the following way. We assume that we can 
prepare the system in the initial state: 

KU = e^^). (17) 

The dynamics of the system is described, in the interac- 
tion picture, by the Hamiltonian: 

H mt = 5J Z1 (18) 
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where 8 is now the constant parameter. Thus, the state 
evolves according to: 



(19) 



and as in the one mode case, after a time T = Air/5 the 
state completes a closed loop. The non-adiabatic phase 
is X(n,n')- Adding decoherence to the system we obtain 
the following phase 



X(n.n') = 4tt(1 - cosa)- 



(20) 
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A 



1 - er~z 
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and at first order in X/S we recover the expression (|11H . 
with Q — 47r(cosa), i.e. the solid angle spanned on the 
parameter sphere in the case of no decoherence. 

Working towards having a realistic description of geo- 
metric phases we have introduced field decoherence in the 
problem of a two-level system interaction with a quan- 
tized field. We analyzed the one and two mode models 
in the adiabatic and non-adiabatic case. We showed that 
when the geometric phase is generated by an adiabatic 
evolution the first correction due to the decoherence of 
the driving field for the no-jump trajectory is only of sec- 
ond order in the decaying rate of the field A. This result 
reinforces the idea that geometric phases can be robust to 
errors, in agreement with previous works which analyze 
the geometric phase under classical noise sources like 
random fluctuations of a classical driving field. We also 
showed that, for the non-adiabatic evolution this is no 
longer the case, and decoherence effects appear already 
in the first order correction term. This result is also in 
accordance with previous works in whic h, again , different 
classical noise sources were considered |l8l I HI . 

Our results are particularly relevant in the experimen- 
tal realizations of these phases, like the one proposed 
in |l5|. and in their use in the implementation of geo- 
metric quantum computation. Understanding the effects 
of decoherence in the geometric evolution of states is the 
first step in finding schemes resilient to this. We are now 
investigating a robust scheme to field decoherence using 
engineered reservoirs. 
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